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Abstract. Prolog's ability to return multiple answers on backtracking 
provides an elegant mechanism to derive reversible encodings of combina- 
torial objects as Natural Numbers i.e. ranking and unranking functions. 
Starting from a generalization of Ackerman's encoding of Hereditarily Fi- 
nite Sets with Urelements and a novel tupling/untupling operation, we 
derive encodings for Finite Functions and use them as building blocks for 
an executable theory of Hereditarily Finite Functions. The more difficult 
problem of ranking and unranking Hereditarily Finite Permutations is 
then tackled using Lehmer codes and factoradics. 

The paper is organized as a self-contained literate Prolog program avail- 
able at http: //logic . csci .unt . edu/tarau/research/2008/pHF F .zip[ 
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1 Introduction 

This paper is an exploration with logic programming tools of ranking and un- 
ranking problems on finite functions and bijections and their related hereditarily 
finite universes. The practical expressiveness of logic programming languages (in 
particular Prolog) are put at test in the process. The paper is part of a larger 
effort to cover in a declarative programming paradigm, arguably more elegantly, 
some fundamental combinatorial generation algorithms along the lines of [13] . 

The paper is organized as follows: section [5] introduces generic ranking/un- 
ranking functions, section [3] introduces Ackermann's encoding in the more gen- 
eral case when urelements are present. Section 3] introduces new tupling/un- 
tupling operations on natural numbers and uses them for encodings of finite 
functions (section [5]), resulting in encodings for Hereditarily Finite Functions 
(section[B]). Ranking/unranking of permutations and Hereditarily Finite Permu- 
tations as well as Lehmer codes and factoradics are covered in section[71 Sections 
[5] and [H] discuss related work, future work and conclusions. 

We will assume that the underlying Prolog system supports the usual higher 
order function-style predicates call/N, findall/3, maplist/N, sumlist/2 or 
their semantic equivalents and a few well known library predicates, used mostly 



for list processing and arithmetics. Arbitrary length integers are needed for some 
of the larger examples but their absence does not affect the correctness of the 
code within the integer range provided by a given Prolog implementation. Oth- 
erwise, the code in the paper, embedded in a literate programming LaTeX file, 
is self contained and runs under SWI- Prolog. Note also that a few utility predi- 
cates, not needed for following the main ideas of the paper, are left out from the 
narrative and provided in the Appendix. 

2 Generic Unranking and Ranking with Higher Order 
Functions 

We will use, through the paper, a generic multiway tree type distinguishing be- 
tween atoms represented as (arbitrary length) integers and subforests represented 
as Prolog lists. Atoms will be mapped to natural numbers in [0. .Ulimit-1]. 
Assuming that Ulimit is fixed, we denote A the set [0 . .Ulimit-1] . We denote 
Nat the set of natural numbers and T the set of trees of type T with atoms in 
A. 

Definition 1 A ranking function on T is a bijection T — )■ Nat. An unranking 
Junction is a bijection Nat — > T. 

Ranking functions can be traced back to Godel numberings |7'8] associated to 
formulae. However, Godel numberings are typically only injective functions, as 
their use in the proofs of Godel's incompleteness theorems only requires injective 
mappings from well-formed formulae to numbers. Together with their inverse 
unranking functions they are also used in combinatorial and uniform random 
instance generation [18113] algorithms. 

2.1 Unranking 

As an adaptation of the unfold operation [9119] . elements of T will be mapped 
to natural numbers with a generic higher order function uiiraiik_ parameterized 
by the the natural number Ulimit and the transformer hmction F: 

unraiik_ (Ulimit, _,N,R) : -NX) ,N<Ulimit , ! ,Ft^. 
uiiraiik_ (Ulimit ,F,N,R) :-N>4Jlimit, 

NO is N-Ulimit, 

call(F,NO,Ns) , 

maplist (unrank_ (Ulimit ,F) ,Ns,R) . 

A global constant provided by the predicate def aultjilimit, will be used 
through the paper to fix the default range of atoms as well as a default unrank 
function: Note also that we will use a syntactically more convenient DCG nota- 
tion, as def aultjilimit will act as a modifier for functional style predicates, 
composed by chaining their arguments automatically with Prolog's DCG trans- 
formation: 



default_ulimit(0)— >[] . 

uiirank(F) — >def ault_ulimit (Ulimit) ,unraiik_ (Ulimit,F) . 



2.2 Ranking 

Similarly, as an adaptation of fold, generic inverse mappings rank_(Ulimit ,G) 
and rcink from T to Nat are defined as: 

rank_ (Ulimit ,_,N,R) : -integer (N) ,N>^,N<Ulimit , ! ,It=N. 

rank_ (Ulimit, G,Ts,R) : -maplist (rank_ (Ulimit ,G) ,Ts,T) , call(G,T,RO) ,R is ROfUlimit . 
rank(G) — >default_ulimit (Ulimit) ,rank_ (Ulimit ,G) . 

Note that the guard in the second definition simply states correctness constraints 
ensuring that atoms belong to the same set A for rank, and uiirank_. This 
ensures that the following holds: 

Proposition 1 // the transformer function F : Nat —^ [Nat] is a bijection with 
inverse G, such that n > ulimit A F{n) = [uq, ...Ui, ...Uk] ^ Ui < n, then unrank 
is a bijection from Nat to T, with inverse rank and the recursive computations 
of both functions terminate in a finite number of steps. 

Proof: by induction on the structure of Nat and T, using the fact that maplist 
preserves bisections. 

3 Hereditarily Finite Sets and Ackermann's Encoding 

The Universe of Hereditarily Finite Sets is best known as a model of the Zermelo- 
Fraenkel Set theory with the Axiom of Infinity replaced by its negation |32|20j . 
In a Logic Programming framework, it has been used for reasoning with sets, 
set constraints, hypersets and bisimulations [6124] . 

The Universe of Hereditarily Finite Sets is built from the empty set (or a set 
of Urelements) by successively applying powerset and set union operations. 

Ackermann's encoding [2|l|llj is a bijection that maps Hereditarily Finite 
Sets {HFS) to Natural Numbers {Nat) as follows: 

f{x) = lix = {} then else Y.aex 2-^^"^ 

Assuming HFS extended with Urelements (atomic objects not having any 
elements) our generic tree representation can be used for Hereditarily Finite 
Sets. 

Ackermann's encoding can be seen as the recursive application of a bijection 
set2nat from finite subsets of Nat to Nat, that associates to a set of (distinct!) 
natural numbers a (unique!) natural number. 



set2nat(Xs,N) : -set2nat (Xs ,0,H) . 



set2nat( [] ,R,R) . 

set2nat([X|Xs] ,Rl,Rn) :-R2 is R1+(1«X) , set2nat (Xs ,R2 ,Rn) . 

With this representation, Ackermann's encoding from HFS to Nat hf s2nat can 
be expressed in terms of our generic rank function as: 

hfs2nat — >def ault_ulimit (Ulimit) ,hf s2nat_ (Ulimit) . 
hf s2nat_ (Ulimit) — >rank_ (Ulimit , set2nat ) . 

where the constant provided by def ault_ulimit controls the segment [0 . . Ulimit-1] 
of Nat to be mapped to urelements. The default value defines "pure" sets, all 
built from the empty set only. 

To obtain the inverse of the Ackerman encoding, we first define the inverse 
nat2set of the bijection set2nat. It decomposes a natural number N into a list 
of exponents of 2 (seen as bit positions equaling 1 in A^'s bitstring representation, 
in increasing order). 

nat2set(N,Xs) : -nat2elements (N,Xs ,0) . 

nat2elements (0 , □ , _K) . 
nat2elements(N,NewEs,Kl) :-N>0, 

B is /\(N,1),N1 is N»1,K2 is Kl+1 ,add_el (B ,K1 ,Es ,NewEs) , 

nat2elements (Nl ,Es ,K2) . 

add_el(0,_,Es,Es) . 
add_el(l,K,Es, [K | Es] ) . 

The inverse of the Ackermann encoding, with urelements in [0 . .Ulimit-1] and 
Ulimit mapped to [] follows: 

nat2hf s_ (Ulimit) — >unrank_ (Ulimit ,nat2set) . 

nat2hfs — >default_ulimit (Ulimit) ,nat2hfs_ (Ulimit) . 

Using an equivalent functional notation, the following proposition summa- 
rizes the results in this subsection: 

Proposition 2 Given id — Xx.x, the following function equivalences hold: 

nat2set o set2nat = id = set2nat o nat2set (1) 

nat2hfs o hf s2nat = id = hfs2nat o nat2hfs (2) 

4 Pairing Functions and Tuple Encodings 

Pairings are bijective functions Nat x Nat Nat. We refer to [S] for a typical 
use in the foundations of mathematics and to [30] for an extensive study of 
various pairing functions and their computational properties. 



4.1 The Pepis-Kalmar- Robinson Pairing Function 

The predicates pepis_pair/3 and pepis_unpair/3 are derived from the function 
pepisJ and its left and right unpairing companions pepis_K and pepis_L that 
have been used, by Pepis, Kalmar and Robinson in some fundamental work on 
recursion theory, decidability and Hilbert's Tenth Problem in |23I10I28] : 

pepis_pair (X,Y,Z) :-pepis_J(X,Y,Z) . 

pepis_unpair (Z,X,Y) :-pepis_K(Z,X) ,pepis_L(Z,Y) . 

pepis_J(X,Y, Z):-Z is ( (1«X) * ( (Y«l)+1) ) -1 . 
pepis_K(Z, X):-Z1 is Z+1 , two_s (Zl ,X) . 

pepis_L(Z, Y):-Z1 is Z+1 ,no_two_s (Zl ,N) , Y is (N-l)»l. 

two_s(N,R) :-even(N) , ! ,H is N»l,two_s(H,T) ,R is T+1 . 
two_s (_ ,0) . 

no_two_s(N,R) :-two_s(N,T) ,R is N // (1«T) . 

even(X) :- =:= /\(1,X) . 

This pairing function given by the formula 

/(x,y) = 2^*(2*2; + l)~l (3) 

is asymmetrically growing, faster on the first argument. It works as follows: 

?- pepis_pair(l , 10,R) . 
R = 41. 

?- pepis_pair(10, 1 ,R) . 
R = 3071. 

?- f indall(R, (between(0,3,A) ,between(0,3,B) ,pepis_pair(A,B,R) ) ,Rs) . 
Rs=[0, 2, 4, 6, 1, 5, 9, 13, 3, 11, 19, 27, 7, 23, 39, 55] 

4.2 Tuple Encodings 

We will now generalize pairing functions to /c-tuples and then we will derive an 
encoding for finite functions. 

The function to_tuple : Nat Nat'' converts a natural number to a fc-tuple 
by splitting its bit representation into k groups, from which the k members in 
the tuple are finally rebuilt. This operation can be seen as a transposition of a 
bit matrix obtained by expanding the number in base 2*^: 

to_tuple(K,N, Ns):- 

Base is 1«K, to_base (Base ,N,Ds) ,maplist (to_maxbits (K) ,Ds ,Bss) , 
mtranspose (Bss ,Xss) , 
maplist (f rom_rbits ,Xss ,Ns) . 



To convert a fc-tuple back to a natural number we will merge their bits, fc at a 
time. This operation uses the transposition of a bit matrix obtained from the 
tuple, seen as a number in base 2*^, with help from bit crunching functions given 
in Appendix: 

f rom_tuple (Nss ,R) :- 

max_bit count (Nss ,L) , length (Nss ,K) .maplist (to_maxbits (L) ,Nss ,Mss) , 
mtranspose (Mss ,Tss) , 

maplist (from_rbits ,Tss ,Ts) , Base is 1«K, f rom_base (Base ,Ts ,R) . 

The following example shows the decoding of 42, its decomposition in bits (right 
to left), the formation of a 3-tuple and the encoding back to 42. 

?- to_tuple(3,42,T) ,to_rbits(2,Bs2) ,to_rbits(l,Bsl) ,from_tuple(T,N) . 
T = [2, 1, 2] , 
Bs2 = [0, 1] , 
Bsl = [1] , 
N = 42 

Note that one can now define pairing functions as instances of the tupling func- 
tions: 

to_pair (N,A,B) :-to_tuple(2,N, [A,B] ) . 
from_pair(X,Y,Z) : -f rom_tuple ( [X, Y] ,Z) . 

One can observe that to_pair and from_pair are the same as the functions 
defined in Steven Pigeon's PhD thesis on Data Compression [25], page 114). 

5 Encoding Finite Functions 

As finite sets can be put in a bijcction with an initial segment of Nat, we can 
narrow down the concept of finite function as follows: 

Definition 2 A finite function is a function defined from an initial segment 
of Nat to Nat. 

This definition implies that a finite function can be seen as an array or a list 
of natural numbers except that we do not limit the size of the representation of 
its values. 

5.1 Encoding Finite Functions as Tuples 

We can now encode and decode a finite function from [0..isr — 1] to Nat (seen as 
the list of its values), as a natural number: 

ftuple2nat ( [] ,0) . 

f tuple2nat (Ns , N) : -Ns=[_ | _] , 

length (Ns,K) ,K1 is K-1, 

f rom_tuple (Ns ,T) ,pepis_pair (Kl ,T, N) . 



nat2ftuple(0, []) . 
nat2f tuple (N, Ns) :-N>0, 

pepis.impair (N,K,F) ,K1 is K+1, 

to_tuple(Kl,F,Ns) . 

As the length of the tuple. K, is usually smaller than the number obtained by 
merging the bits of the K-tuple, we have picked the Pepis pairing function, ex- 
ponential in its first argument and linear in its second, to embed the length of 
the tuple needed for the decoding. The encoding/decoding works as follows: 

?- ftuple2nat( [1,0,2,1,3] ,M). 
N = 21295 

?- nat2f tuple (21295, T) . 
T=[l,0,2,l,3] 

?-ints_f rom(0 , 15, Is) ,maplist (nat2f tuple , Is ,Ts) . 
Ts=[[0] , [0,0] , [1] , [0,0,0] , [2] , [1,0] , [3] , 

[0,0,0,0] , [4] , [0,1] , [5] , [1,0,0] , [6] , 

[1,1] , [7] , [0,0,0,0,0]] 

Note that 
nat (0) . 

nat(N) :-nat(Nl) ,N is Nl+1. 

iterative_fun_generator (F) :-nat(N) ,nat2f tuple (N,F) . 
provides an iterative generator for the stream of finite functions. 

5.2 Deriving Encodings of Finite Functions from Ackermann's 
Encoding 

Given that a finite set with n elements can be put in a bijection with [0..iV— 1], a 
finite functions / : [0..n— 1] Nat can bo represented as the list [/(0).../(n— 1)]. 
Such a list has however repeated elements. So how can we turn it into a set with 
distinct elements, bijectively? 

The following two predicates provide the answer. 

First, wc just sum up the list of the values of the function, resiilting in a 
monotonically growing sequence (provided that we first increment every number 
by 1 to ensure that values do not break monotonicity). 

f un2set ([],[]). 

f un2set ( [A | As] , Xs ) : -f indall (X , pref ix.sum ( A , As , X) , Xs ) . 

pref ix_sum(A, As ,R) : -append (Ps ,_, As) , length (Ps,L) , 
sumlist(Ps,S) ,R is A+S+L. 

The inverse of f un2set reverting back from a set of distinct values collects the 
increments from a term to the next (and ignores the last one): 



set2fmi([] , []) . 

set2f un ( [X | Xs] , [X | Fs] ) : -set2f un (Xs , X , Fs) . 
set2fun( [],_,[]). 

set2fun([X|Xs] ,Y, [A|As]) :-A is (X-Y)-l , set2fun(Xs ,X , As) . 

Proposition 3 The following Junction equivalences hold: 

fun2set o set2fun = id = set2fun o fun2set (4) 
The mapping and its inverse work as follows: 

?- fun2set( [1 ,0,2, 1,2] ,Set) , set2f un(Set ,Fun) . 
Set = [1, 2, 5, 7, 10] , 
Fun = [1, 0, 2, 1, 2] . 

By combining this bijection with Ackermann's encoding's basic step set2nat 
and its inverse nat2set, we obtain an encoding from finite functions to Nat as 
follows (with DCG notation used to express function composition): 

nat2fun — > nat2set , set2fuii. 
f un2nat — > f uii2set , set2nat . 

?- nat2fun(2008,F) ,fun2nat(F,N) . 
F= [3, 0, 1, 0, 0, 0, 0], N= 2008 

Proposition 4 The following function equivalences hold: 

nat2fun o fun2nat = id = fun2nat o nat2fun (5) 

One can see that this encoding ignores Os in the binary representation of a 
number, while counting 1 sequences as increments. Alternatively, Run Length 
Encoding of binary sequences [21j encodes Os and Is symmetrically, by counting 
the numbers of Is and Os. This encoding is reversible, given that Is and Os 
alternate, and that the most significant digit is always 1: 

bits2rle([] ,[]):-!. 
bits2rle([_] , [0]) :-! . 

bits2rle([X,Y|Xs] ,Rs) :-X=Y, ! ,bits2rle ( [Y | Xs] , [C|Cs]) ,C1 is C+l , Rs=[Cl | Cs] . 
bits2rle ( [_ | Xs] , [0 | Rs] ) : -bits2rle (Xs , Rs) . 

rle2bits([] , [] ) . 

rle2bits([N|Ns] ,NewBs) : -rle2bits (Ns ,Xs) , 
( []=Xs->B is 1 
; Xs=[Xl I _] ,B is 1-Xl 
), 

Nl is N-|-l,ndup(Nl,B,Bs) ,append(Bs,Xs,NewBs) . 



By composing bits2rle and rle2bits with converters to/from bitlists, we ob- 
tain the bijcction nat2rle : Nat ~> [Nat] and its inverse rle2nat : [Nat] —^ Nat 

nat2rle — > to_rbitsO ,bits2rle . 
rle2nat — > rle2bits ,from_rbits . 

to_rbitsO(0, [] ) . 

to_rbitsO(N,R) : -N>0, to_rbits (N,R) . 

Proposition 5 The following function equivalences hold: 



6 Encodings for "Hereditarily Finite Functions" 

One can now build a theory of "Hereditarily Finite Functions" [HFF) cen- 
tered around using a transformer like nat2f tuple, nat2fun, nat2rle and 
ftuple2nat, fun2nat, rle2nat in way similar to the use of nat2set and 
set2nat for HFS, where the empty function (denoted [] ) replaces the empty set 
as the quintessential "urfunction" . Similarly to Urelements in the HFS theory, 
"urfunctions" (considered here as atomic values) can be introduced as constant 
functions parameterized to belong to [O..Ulimit — 1]. 



By using the generic unrank_ and rank predicates defined in section[2]we can 
extend the bijections defined in this section to encodings of Hereditarily Finite 
Functions. By instantiating the transformer function in unrank_ to nat2f tuple, 
nat2f un and nat2rle we obtain (with DCG notation expressing composition of 
functional predicates): 

nat2hff — > def ault_ulimit (D) ,nat2hf f _ (D) . 
nat2hffl — > def ault_ulimit (D) ,nat2hf f 1_ (D) . 
nat2hf f 2 — > def ault_ulimit (D) ,nat2hf f 2_ (D) . 

nat2hf f _ (Ulimit) — > unrank_ (Ulimit ,nat2fun) . 
nat2hff l_(Ulimit) — > unrank. (Ulimit ,nat2f tuple) . 
nat2hff2_ (Ulimit) — > unrank_ (Ulimit ,nat2rle) . 

By instantiating the transformer function in rank we obtain: 

hff2nat — > rank(fun2nat) . 
hff2natl — > rank(ftuple2nat) . 
hff2nat2 — > raiik(rle2nat) . 

The following examples show that nat2hff, nat2hff 1 and nat2hff2 are 
indeed bijections, and that the resulting HFF-tTees are typically more compact 
than the HFS-tiee associated to the same natural number. 

?- nat2hff (42,H) ,hff2nat(H,N) . 



nat2rle o rle2nat = id = rle2nat o nat2rle 



(6) 



H = [[[]] , [[]] , [[]]] , 



N = 42 



?- nat2hff 1(42,H) ,hff2natl(H,N) . 
H = [[[[] , [] , []] , []]] , 
N = 42 

?- nat2hff2(42,H) ,hff2nat2(H,N) . 
H = [[] , [] , [] , [] , [] , []] , 
N = 42 

Note that 

?-iiat(N) ,iiat21iff (N,HFF) . 
?-iiat(N) ,nat2hff 1(N,HFF) . 
?-iiat(N) ,iiat2hff2(N,HFF) . 

provide iterative generators for the (recursively enumerable!) stream of heredi- 
tarily finite functions. 

The resulting HFF with urfunctions (seen as digits) can also be used as gen- 
eralized numeral systems with possible applications to building arbitrary length 
integer implementations. 

?- nat2hff_(10,1234567890,HFF) . 
[3, 2, 0, 1, 7, 0, 1, 2, 0, 2, 2] 

Proposition 6 The following function equivalences hold: 

nat2hffl o hff2natl = id = hff2natl o nat2hffl (7) 

nat2hf f o hf f2nat = id = hf f2nat o nat2hf f (8) 

7 Encoding Finite Bijections 

To obtain an encoding for finite bijections (permutations) we will first review a 
ranking/unranking mechanism for permutations that involves an unconventional 
numeric representation, factoradics. 

7.1 The Factoradic Numeral System 

The factoradic numeral system |14j replaces digits multiplied by power of a 
base N with digits that multiply successive values of the factorial of N. In the 
increasing order variant f r the first digit do is 0, the second is di S {0, 1} and 
the A^-th is djv G [0..A'' — 1]. The left-to-right, decreasing order variant fl is 
obtained by reversing the digits of f r. 



?- fr(42,R) ,rf (R,N) . 
R = [0, 0, 0, 3, 1] , 
N = 42 



?- fl(42,R) ,lf (R,N) . 
R = [1, 3, 0, 0, 0] , 
N = 42 

The Prolog predicate fr handles the special case for and calls frl which 
recurses and divides with increasing values of N while collecting digits with mod: 

factoradics of N, right to left 
fr(0, [0]) . 

fr(N,R) :-N>0,frl(l,N,R) . 
frl(_,0,[]). 

frl(J,K, [KMJ|Rs]) :-K>0,KMJ is K mod J,J1 is J+1,KDJ is K // J, 
frl(Jl,KDJ,Rs) . 

The reverse f 1, is obtained as follows: 
fl(N,Ds) :-fr(N,Rs) .reverse (Rs ,Ds) . 

The predicate If (inverse of f l) converts back to decimals by summing up results 
while computing the factorial progressively: 

If (Ls,S) :-length(Ls,K) ,K1 is K-1 , If (Kl ,_ ,S ,Ls , [] ) . 

7, from list of digits of factoradics, back to decimals 
If (0,1,0)— >[0] . 

If (K,N,S)— >[D] ,{K>0,K1 is K-l> , If (Kl ,N1 ,S1) , {N is K*N1,S is Sl-fD*N}-. 
Finally, rf , the inverse of f r is obtained by reversing f 1. 
rf (Ls,S) : -reverse (Ls.Rs) ,lf (Rs,S) . 



7.2 Ranking and unranking permutations of given size with Lehmer 
codes and factoradics 

The Lehmer code of a permutation / is defined as the number of indices j such 
that 1 < j < i and /(j) < f{i) [T^. 

Proposition 7 The Lehmer code of a permutation determines the permutation 
uniquely. 

The predicate perm2nth computes a rank for a permutation Ps of Size>0. It 
starts by first computing its Lehmer code Ls with perni_lehmer. Then it asso- 
ciates a unique natural number N to Ls, by converting it with the predicate If 
from factoradics to decimals. Note that the Lehmer code Ls is used as the list 
of digits in the factoradic representation. 



perm2nth (Ps , Size , N) : - 

length(Ps , Size) .Last is Size-1, 
ints_from(0,Last ,Is) , 
perm_lehmer(Is,Ps,Ls) , 
If (Ls,N) . 

The generation of the Lehmer code is surprisingly simple and elegant in Prolog. 
We just instrument the usual backtracking predicate generating a permutation to 
remember the choices it makes, in the auxiliary predicate select _and_remember! 

associates Lehmer code to a permutation 
perm_lehmer ([],[],[]). 
penn_lehmer (Xs , [X | Zs] , [K | Ks] ) : - 

select _and_remember (X,Xs ,Ys ,0,K) , 
perm_lehmer (Ys , Zs , Ks) . 

y. remembers selections - for Lehmer code 
select_aiid_remember (X, [X|Xs] ,Xs,K,K) . 

select_and_remember(X, [Y|Xs] , [Y| Ys] ,K1,K3) :-K2 is Kl+1, 
select_and_remember (X , Xs , Ys , K2 , K3) . 

The predicate nat2perm provides the matching unranking operation associ- 
ating a permutation Ps to a given Size>0 and a natural number N. 

nth2perm(Size,N, Ps):- 
fl(N,Ls) , length (Ls,L) , 

K is Size-L,Last is Size-1 , ints_from(0, Last , Is) , 

zeros(K,Zs) ,append(Zs,Ls,LehmerCode) , 
perm_lehmer (Is.Ps.LehmerCode) . 

Note also that perm_lehmer is used (reversibly!) this time to reconstruct the 
permutation Ps from its Lehmer code. The Lehmer code is computed from the 
permutation's factoradic representation obtained by converting N to Ls and then 
padding it with O's. One can try out this bijcctive mapping as follows: 

?- nth2perin(5,42,Ps) ,perin2nth(Ps, Length, Nth) . 
Ps = [1, 4, 0, 2, 3] , 
Length = 5, 
Nth = 42 

?- nth2perm(8,2008,Ps) ,perin2nth(Ps, Length, Nth) . 
Ps = [0, 3, 6, 5, 4, 7, 1, 2] , 
Length = 8, 
Nth = 2008 

7.3 A bijective mapping from permutations to Nat 

One more step is needed to to extend the mapping between permutations of 
a given length to a bijcctive mapping from/to Nat: we will have to "shift to- 
wards infinity" the starting point of each new bloc of permutations in Nat as 
permutations of larger and larger sizes are enumerated. 



First, we need to know by how much - so we compute the sum of all factorials 
up to iV!. 

'/, fast computation of the sum of all factorials up to N! 
sf (0,0) . 

sf (N,R1) :-N>0,Nl is N-1 ,ndup(Nl , 1 ,Ds) ,rf ( [0 | Ds] ,R) ,R1 is R+1 . 

This is done by noticing that the factoradic representation of [0,1,1,..] does just 
that. The stream of all such sums can now be generated as usual: 

sf (S) :-nat(N) ,sf (N,S) . 

What we are really interested into, is decomposing N into the distance to the 
last sum of factorials smaller than N, N_M and its index in the sum, K. 

to_sf(N, K,N_M) :-nat(X) ,sf (X,S) ,S>N, ! ,K is X-1 , sf (K ,M) ,N_M is N-M. 

Unranking of an arbitrary permutation is now easy - the index K determines the 
size of the permutation and N_M determines the rank. Together they select the 
right permutation with nth2perm. 

nat2perm(0, [] ) . 

nat2perm(N,Ps) :-to_sf (N, K,N_M) ,nth2perm(K,N_M,Ps) . 

Ranking of a permutation is even easier: we first compute its Size and its rank 
Nth, then we shift the rank by the sum of all factorials up to Size, enumerating 
the ranks previously assigned. 

perm2nat ( [] ,0) . 

perm2nat(Ps,N) : -perm2nth (Ps , Size, Nth) ,sf (Size, S) ,N is S-RJth. 

?- nat2perm(2008,Ps) ,perm2nat (Ps ,N) . 
Ps = [1, 4, 3, 2, 0, 5, 6] , 
N = 2008 

As finite bijections are faithfully represented by permutations, this construction 
provides a bijection from Nat to the set of Finite Bijections. 

Proposition 8 The following function equivalences hold: 

nat2perm o perm2nat = id = perm2nat o nat2perm (9) 

7.4 Hereditarily Finite Permutations 

By using the generic unrank_ and rank predicates defined in section [5] we can 
extend the nat2perm and perm2nat to encodings of Hereditarily Finite Permu- 
tations (HFP). 

nat2hfp — > def ault_ulimit (D) ,nat2hf p_ (D) . 
nat2hf p_ (Ulimit) — > unrank_ (Ulimit ,nat2perm) . 
hfp2nat — > raiik(perm2nat) . 



The encoding works as follows: 



?- nat2hfp(42,H) ,hfp2nat(H,N) .write (H) ,nl. 

H= [[], [[], [[]]], [[[]], []], [[]], [[], [[]], [[], [[]]]]], 
N = 42 

Proposition 9 The following function equivalences hold: 

nat2hfp o hfp2nat = id = hfp2nat o nat2hfp (10) 

8 Related work 

Natural Number encodings of Hereditarily Finite Sets have triggered the interest 
of researchers in fields ranging from Axiomatic Set Theory and Foundations of 
Logic to Complexity Theory and Combinatorics |32lllll2lll4l20ll6l31l3j . Com- 
putational and Data Representation aspects of Finite Set Theory have been 
described in logic programming and theorem proving contexts in [6124122] . Pair- 
ing functions have been used work on decision problems as early as [23110127129] . 
The tuple functions we have used to encode finite functions are new. While fi- 
nite functions have been used extensively in various branches of mathematics 
and computer science, we have not seen any formalization of hereditarily Finite 
Functions or Hereditarily Finite Bijections as such in the literature. 

9 Conclusion and Future Work 

We have shown the expressiveness of logic programming as a metalanguage for 
executable mathematics, by describing natural number encodings, tupling/untu- 
pling and ranking/unranking functions for finite sets, functions and permuations 
and by extending them in a generic way to Hereditarily Finite Sets, Hereditarily 
Finite Functions and Hereditarily Finite Permutations. 

In a Genetic Programming context [15126] . the bijections between bitvec- 
tors/natural numbers on one side, and trees/graphs representing HFSs, HFFs, 
HPPs on the other side, suggest exploring the mapping and its action on various 
transformations as a phenotype-genotypc connection. 

We also foresee interesting applications in cryptography and steganography. 
For instance, in the case of the permutation related encodings - something as 
simple as the order of the cities visited or the order of names on a greetings 
card, seen as a permutation with respect to their alphabetic order, can provide 
a steganographic encoding/decoding of a secret message by using predicates like 
nat2perm and perm2nat. 

Last but not least, the use of a logic programming language to express in a 
generic way some fairly intricate combinatorial algorithms predicts an interesting 
new application area. 
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A Appendix 

To make the code in the paper fully self contained, we list here some auxiliary 
functions. 

Integer list operations These are some simple utility predicates: 
7. generates integers From. .To 

ints_from(From,To,Is) :-f in<iall(I,between(From,To,I) ,Is) . 

y. replicates X, N times 
ndupCO, _,[]). 

ndupCN.X, [X|Xs]) :-M>0,Nl is N-l,ndup(Nl,X,Xs) . 
zeros (N,Zs) :-ndup(N,0,Zs) . 

Matrix Transposition This code transposes a matrix represented as list of lists, 
mtranspose ([],[]):-!. 

mtranspose ( [Xs] , Css) : - ! , to_coliimns (Xs , Css) . 
mtranspose ( [Xs | Xss] , Css2) : - ! , 

mtranspose (Xss ,Cssl) , 

to_colnmns (Xs , Cssl , Css2) . 



to_col\imns ( [] , [] ) . 



to_columns ( [X | Xs] , [ [X] | Zs] ) : -to_columns (Xs , Zs) . 
to_coliimns ( [] , Css , Css) . 

to.columns ( [X I Xs] , [Cs|Cssl] , [[X|Cs] |Css2]) :- to_columiis(Xs,Cssl,Css2) . 

Bit crunching functions The following functions implement conversion oper- 
ations between bitlists and numbers. Note that our bitlists represent binary 
numbers by selecting exponents of 2 in increasing order (i.e. "right to left"). 

7. conversion to list of digits in given base 
to_base (Base,N,Bs) :-to_base(N,Base,0,Bs) . 

to_base(N,R,_K,Bs) :-N<R,Bs=[N] . 
to_base(N,R,K, [B|Bs]) :-N>*, 

B is N mod R, Nl is N//R,K1 is K+1, 

to.base (Nl,R,Kl,Bs) . 

7. conversion from list of digits in given base 
f rom_base (_Base , [] , 0) . 

f rom_base (Base, [X I Xs] ,N) :-f roin_base (Base ,Xs ,R) ,N is X-|-R*Base. 

■/, conversion to list of bits, right to left 
to_rbits(N,Bs) :-to_base(2,N,Bs) . 

7. conversion from list of bits, right to left 
f rom_rbits (Bs , N) : -f rom_base (2 , Bs , N) . 

7o counting how many bits a niimber needs 
bitcovmt (N , K) : -Ifc=<l , K=l . 

bitcount(N,K) :-M>l,Nl is N»l ,bitcount (Nl ,K1) ,K is Kl+1. 
'/. finds the larges bitcount for a list 

max_bitcoimt(Nss,L) :-maplist (bitcount, Nss,Ls) ,mcix_list(Ls,L) . 

7. pads up to maxbits , if needed 
to_maxbits (Maxbits ,N,Rs) :- 

to_base(2,N,Bs) ,length(Bs,L) ,ML is Maxbits-L, 

ndup (ML , , Zs) , append (Bs , Zs , Rs) . 



